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ABSTRACT 



The Simultaneous Medicina- P/anc/: Experiment (SiMPlE) is aimed at observing 
a selected sample of 263 extragalactic and Galactic sources with the Medicina 32-m 
single dish radio telescope in the same epoch as the Planck satellite observations. 

The data acquired with a frequency coverage down to 5 GHz, also combined with 
Planck at frequencies above 30 GHz, will constitute a useful reference catalogue of 
\ bright sources over the whole Northern hemisphere. Furthermore, source observations 

■ performed in different epochs and comparison with other catalogues allow the inves- 

tigation of source variabilities on different timescales. 

In this work, we describe the sample selection, the on-going data acquisition cam- 
paign, the data reduction procedures, the developed tools, and the comparison with 
other data-sets. 

We present the data at 5 and 8.3 GHz for the SiMPlE Northern Sample consisting 
of 79 sources with 5 ^ 45° selected in our catalogue and observed during the first 
6 months of the project. A first analysis of their spectral behaviour and long-term 
variability is also presented. 



Key words: galaxies: active - radio continuum: galaxies - radio continuum: general 
- cosmic microwave background. 



1 INTRODUCTION 

Microwave band observations with satellites mainly aimed 
at the study of the cosmic microwave background (CMB) 
also offer a unique opportunity to investigate the spectral 
properties of radio sources in a poorly explored frequency 
range, partially inaccessible from ground. 

The study of the Spectral Energy Distribution (SED) 
and variability of sources from radio to far-infrared (FIR) 
wavelengths is crucial to probe the physics of the innermost 
regions of Active Galactic Nuclei (AGN) and radio galax- 
ies and of some classes of variable Galactic radio sources 
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(i.e. microquasars and binary active systems). These objects 
show remarkable outbursts and studying their SED simulta- 
neously from radio to FIR in several epochs helps constrain- 
ing the outburst mechanisms and the processes that activate 
their radio emission. 

The project described in this paper aims at accumulat- 
ing multifrequency radio observations of the northern sky 
at 5, 8.3 and 22 GHz almost simultaneous with those taken 
by the ESA Planck satellite (Tauber et al. 2010a) from 30 
to 857 GHz. This will allow the construction of the SED 
of a significant sample of sources on a wide frequency range 
through almost coeval data. The simultaneity of the observa- 
tions is crucial to evaluate the contribution from variability 
and spectral effects in affecting the physical interpretation 
of source observational properties. 
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Catalogues based on the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite maps have been 
produced using several source extraction techniques 
(Lopez-Caniego et al. 2006; Chen & Wright 2009; Bennett 
et al. 2003; Hinshaw et al. 2007; Gold et al. 2011). Theh 
completeness limit (i.e. the flux density level above which 
all the sources in the observed area are listed in the 
catalogue) is typically ~ IJy at 23 GHz. Massardi et al. 
(2009) have combined blind and non-blind approaches (see 
Lopez-Caniego et al. 2007, Gonzalez-Nuevo et al. 2008) 
exploiting the Mexican Hat Wavelets 2 (MHW2) fiher 
(Gonzalez-Nuevo et al. 2006) to extract an all-sky catalogue 
of 516 sources with |6| > 5°. It is almost 91 per cent 
complete above 1 Jy and constitutes the New Extragalactic 
WMAP Point Sources (NEWPS) catalogue. Almost all 
(484) sources in the sample were previously catalogued as 
extragalactic (457) or Galactic (27) objects. The remaining 
32 candidate sources do not have counterparts in lower- 
frequency all-sky surveys with comparable flux densities 
and may therefore be just high peaks in the distribution 
of other components present in the maps. If they are all 
spurious, the reliability of the sample (i.e. the probability 
that a source listed in the sample is a genuine radio source) 
is 93.8 per cent. 

To date, the Australia Telescope 20-GHz Survey 
(AT20G, Murphy et al. 2010, Massardi et al. 2010) pro- 
vides the deepest complete ground-based sample of the high- 
frequency Southern sky. It is a blind survey performed with 
the Australia Telescope Compact Array (ATCA) in the years 
2004 — 2008. Its final catalogue consists of 5890 sources 
above a flux limit of 40 mjy and is an order of magnitude 
larger than any previous catalogue of high-frequency radio 
sources. No analogous sample, deeper than the NEWPS cat- 
alogue has been observed at frequencies above 8.4 GHz in 
the Northern hemisphere. 

The ESA Planck satellite (Tauber et al. 2010a) is pro- 
viding a blind survey of the entire sky in nine frequency 
bands (30, 44, 70, 100, 143, 217, 353, 545, 857 GHz), with 
Full Width Half Maximum (FWHM) resolution ranging 
from 33 to 5 arcmin (Mennella et al. 2011, Planck HFI Core 
Team et al. 2011) and at several epochs. Leach et al. (2008) 
estimated the possibilities of the MHW2-filter-based detec- 
tion techniques applied to the Planck satellite maps, consid- 
ering two all-sky surveys. They found that the expected de- 
tection limits range from ~ 0.4 Jy at 30 GHz to ~ 0.22 Jy at 
100 GHz, which implies that all the NEWPS sources should 
be detectable also in the Planck maps up to ~ 100 GHz. 

The Planck Early Release Compact Source Catalogue 
(ERCSC) has been recently released (The Planck Collabo- 
ration 2011a) on the basis of more than one full coverage of 
the entire sky. The catalogue refers, in fact, to a stage when a 
first full-sky map plus a second run on 60 per cent of the sky 
were available. A Monte-Carlo algorithm was implemented 
to select reliable sources among all the extracted candidates 
such that the cumulative reliability of the catalogue is 90 
per cent. There is no requirement on completeness for the 
ERCSC. The lOo" photometric flux density limit of the cat- 
alogue at |6| > 30° is 0.49, 1.0, 0.67, 0.5, 0.33, 0.28, 0.25, 
0.47 and 0.82 Jy at each of the nine frequencies between 30 
and 857 GHz. 

The Planck ERCSC "provides a robust list of stars with 
dust shells, stellar cores, radio galaxies, blazars, infrared 



luminous galaxies. Galactic interstellar medium features, 
915 cold molecular cloud core candidates, 189 Sunyaev- 
Zel'dovich cluster candidates as well as unclassifled sources" 
(The Planck Collaboration 2011a). Many of them are ob- 
ject of dedicated papers about Galactic science, extragalac- 
tic sources, and Sunyaev-Zel'dovich effects and cluster prop- 
ertied. 

The source list, with more than 15000 unique sources, is 
ripe for follow-up characterisation with Herschel and several 
ground-based observing facilities. 

The first version of the Planck Legacy Catalog, together 
with other Planck products and a first set of cosmological 
papers, will be released in the first months of 2013, i.e. at 
the end of the proprietary period of the data from the first 
two surveys (namely, 15 months of observation). 

Planck observations help identifying the properties of 
the SED in the high frequency radio and FIR bands (The 
PZancfc collaboration 2011b, c), and benefit from the ground- 
based observations at longer wavelengths for the reconstruc- 
tion of the properties of different emission components to 
study the spectral behaviour of different radio source popu- 
lations . 

An analysis of the variability on the five-year WMAP 
point sources shows that a high fraction of the sources are 
variable at more than 99 per cent confidence, and these are 
in general the brighter ones (Wright et al. 2009). In gen- 
eral, most AGNs are variable at these radio frequencies. 
Their long-term variability has been well studied over the 
years at similar or higher frequencies, for example, by the 
University of Michigan Radio Astronomy Observatory (UM- 
RAO) and Metsahovi Radio Observatory research teams for 
sources in the northern sky and equatorial regions (Hov- 
atta et al. 2007, Hughes et al. 1992). Also, the comparison 
between the WMAP catalogue and the Planck ERCSC at 
the corresponding bands, that can be performed in a gen- 
eral, statistical sense because of source variability, reveals no 
systematic difference between the WMAP and ERCSC flux 
densities, while the signiflcant scatter confirms that variabil- 
ity is an issue (The Planck collaboration 2011a). Variability 
enhances the scattering of flux density comparisons. It also 
results in a bias in favour of sources in a bright phase at the 
selection epoch. This will make some of the variable (and 
mostly bright flat spectrum) ones appear with a Gigahertz 
Peaked-Spectrum (GPS) like distorted spectrum character- 
ising flaring phases. 

Sadler et al. (2006) estimated a median debiased vari- 
ability of 6.9 per cent at 20 GHz and on timescales of a year 
for a 100-mJy flux density limited sample of extragalactic 
radio sources, with only a few sources more variable than 30 
per cent. Most recently, Massardi et al. (2011) estimated at 
18 GHz a median debiased variability on 9 months of 9 per 
cent for a sample of AT20G sources with ^oghz >500 mJy. 

The flux density levels achieved during outbursts by the 
above mentioned classes of Galactic radio sources are high 
enough to be observed with good signal-to-noise ratio (SNR) 
by Planck. 

Driven by all the above reasons, several projects are 



The Planck Early Results papers are pub- 
licly available at the following web address 
jhttp: //www. sciops ■ esa. Int/index ■php?project=PLANCKfepage=Planck_Publlshed. 
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Figure 1. The SiMPlE sample: the Northern Sample (diamonds) 
covers the area with 6 ^ 45°; asterisks show the Southern Sample; 
squares indicate the central positions of the mapped sky patches. 
The dotted lines represent a 10-degrees masked band centred on 
the Galactic plane (middle dotted line). 



carrying out observations simultaneous with the Planck 
satellite at various frequency bands. They involve observa- 
tional facilities such the Australia Telescope Compact Ar- 
ray (ATCA), the Effelsberg Radiotelescope (Germany), the 
Institut de Radioastronomie Millimetrique (IRAM) 30 m 
Telescope (Spain) the Metsahovi Radio Observatory (Fin- 
land), the Very Large Array / Expanded Very Large Array 
(VLA/EVLA) (USA), operating in the frequency range be- 
tween 4.5 and 40 GHz, 2.6 and 43 GHz, and at 86.2, 142.3 
GHz, at 37 GHz, and at 5, 8, 22, 43 GHz, respectively (see 
The Planck CoUaboration 2011b, c). 

The Planck ATCA Coeval Observations (PAGO, Mas- 
sardi et al. 2011) project is the largest one and has followed 
up a sample of 482 AT20G extragalactic sources in the fre- 
quency range between 4.5 and 40 GHz in the period between 
July 2009 and August 2010. Several sources were observed 
more than once to study variability. 

We present here the Simultaneous Medicina Planck 
Experiment (SiMPlE), designed to complement the PAGO 
project in the Northern hemisphere using the Medicina 32-m 
single dish to observe a sample of 263 sources at 5, 8.3 and 
22 GHz almost simultaneously to the Planck observations. 

The SiMPlE data alone at low frequencies and com- 
bined with Planck for frequencies above 30 GHz, constitute 
a useful reference catalogue of bright sources over the whole 
Northern hemisphere. The sample selection criteria are de- 
scribed in Section [2l the observing strategy is given in Sec- 
tion [31 and data reduction techniques are discussed in Sec- 
tion|4l Data and spectral analysis for the 79 NEWPS sources 
with 5 J5 45° at 5 and 8.3 GHz are presented and discussed 
in Section [S] Finally, the description of the present status of 
the project, our main findings, and some future perspectives 
are summarised in Section (6] 



2 SAMPLE SELECTION 

The selection of our sample of sources largely exploits the 
NEWPS sample (Massardi et al. 2009). At the time of plan- 
ning out observations it constituted a unique opportunity 
to select high-frequency samples of bright sources in the 
Northern hemisphere, including 252 sources with declina- 
tion (5 > 0°. 

Cross-matches of the northern NEWPS sources with 
low-frequency wide-area surveys like the GB6 (Gregory et 
al. 1996) catalogue at 4.85 GHz or the NVSS (Condon et 
al. 1998) at 1.4 GHz have identified the candidate counter- 
parts for 234 sources, 11 of which show multiple associations. 
Identifications based on the NASA Extragalactic Database 
(NED) classified 225 of them as extragalactic objects and 
the remaining 9 as Galactic. 18 sources do not have coun- 
terparts in lower- frequency all-sky surveys with comparable 
flux densities and may be spurious sources. 

Because of the large WMAP beam and, as a conse- 
quence, of the positional uncertainty of the NEWPS detec- 
tions, in planning follow-up observations we have used the 
coordinates of the low-frequency counterpart, where avail- 
able, as target positions. In case of spurious sources or mul- 
tiple associations we have chosen the WMAP positions as 
centre of a region that we have mapped at 5 GHz with the 
Medicina radio telescope in order to identify, in the case it 
was a compact radio source, the low- frequency counterpart 
of the WMAP detection (see Sect. [331 for further details). 
The SiMPlE observations will therefore help to assess pre- 
cisely the reliability of the detection techniques that have 
been applied to the WMAP maps. 

Comparisons of the NEWPS sample detection with the 
AT20G Bright Source Sample (Massardi et al. 2008) data 
indicate that the NEWPS sample is almost 91 per cent com- 
plete down to 1 Jy. The completeness of the sample is deter- 
mined by the flux density scales at the epoch of the WMAP 
observations, and variability can alter it. However, the com- 
parison with the coeval 30-GHz Planck observations, that 
are achieving deeper flux densities, will allow us to quantify 
the level of completeness of our sample at least at 22 GHz, 
once we will properly account for the difference in frequency. 
The 22-GHz data and the comparisons with the Planck re- 
sults will be presented in future papers from our group. 

In addition to the Galactic sources in the NEWPS cata- 
logue at 5 > 0° , the SiMPlE project includes the monitoring 
of a sample of 10 sources representative of various classes 
of highly variable objects, in case they may show a strong 
outburst during the Planck mission. Massive X-ray binary 
systems, in fact, alternate quiescent periods and strong out- 
bursts: about 10 per cent of the more than 200 known binary 
systems are radio- loud (Mirabel & Rodrigues 1999). Some 
of them (e.g. Cyg X-3) are expected to reach flux densities 
up to 20 Jy (Trushkin et al., 2008; Szostek et al., 2008) at 
8.4 GHz. Luminous Blue Variable stars (e.g. Eta Carinae) 
show sudden outbursts during which the flux density can in- 
crease to several Jy at centimetre wavelengths and tens of Jy 
in the millimetre band. Finally, active binary stars (e.g. RS 
CVn and Algol) have periods of quiescence, with flux den- 
sities of a few tens of mJy, and active periods characterised 
by flares, lasting several weeks, every 2-3 months, reaching 
from several hundred mJy up to some Jy at cm wavelengths 
(Umana et al. 1995). Their radio spectra show the maximum 
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of emission at frequencies higher than 10 GHz and, during 
the impulsive phases, this moves up to about 100 GHz. 

Pre-emptive observations of these objects began in 2009 
with the Nolo telescope at 43 GHz and continued in Medic- 
ina, starting from June 2010 in the framework of the SiM- 
PIE experiment at lower frequencies. These observations are 
on-going and will cover at least 1 years of Planck data ac- 
quisition. The results and the comparison with Planck data 
as well as 22-GHz data will be shown in future papers. 

In summary, the SiMPlE sample contains 263 targets 
of which 253 are extracted from the NEWPS sample. 29 of 
these have been observed by mapping the surrounding re- 
gion, while cross- On- The- Fly (OTF) scans have been used 
to observe the others. Fig. [T] shows the distribution on the 
Northern hemisphere of the observed sample. 79 objects 
have declination 45° and constitute the so-called 'North- 
ern sample' for which in Section [S] we present the results of 
the observations in the epochs between June and Novem- 
ber 2010. The remaining epochs for the full sample will be 
discussed in future papers. 



3 OBSERVING STRATEGY WITH THE 
MEDICINA RADIO TELESCOPE 

3.1 Simultaneity with Planck 

The Planck satellite scans the sky in circles passing close to 
the ecliptic poles and covers the whole sky in about 7 months 
(Dupac & Tauber 2005; see also section 1 of Mandolesi et 
al. 2010 and references therein). All the beams (Tauber et al 
2010b) of the entire Planck field of view cross the same sky 
position in a few days at low and intermediate ecliptic lati- 
tudes, in weeks or few months at high ecliptic latitudes, and 
then pass again on it at the subsequent survey, so allowing to 
recover source flux density variations on different timescales 
(Burigana 2000, Terenzi et al. 2002, 2004). The extraction 
of the source flux density at Planck frequency channels is 
based on the analysis of frequency maps produced by op- 
timally weighting the Time Ordered Data (TOD) in pixel 
space. Therefore, we can not achieve information from them 
on the very short term variability, a study that necessarily 
requires to exploit the information contained in time-ordered 
data, less sensitive than the channel maps. We then consider 
observations to be 'simultaneous' with the satellite if per- 
formed within 10 days from the satellite observations at any 
of its frequencies. In practice, this has been also a reasonable 
compromise with the scheduled observational days available 
to the project. This typical time sampling does not prevent 
the extraction of relevant variability information at least for 
the large majority of extragalactic sources, whose flux den- 
sities change appreciably on significantly longer timescales. 
The possibility to carry out observations, possibly simul- 
taneously with Planck, of bright sources (in particular of 
objects expected to show faster variability) with a finer 
timescale at the Medicina radio telescope is under inves- 
tigation for the next SiMPlE campaigns. In general, the 
Planck On-Flight Forecaster (Massardi & Burigana 2010) 
was applied to predict when our target sources were being 



Table 1. Observational details used to plan the cross scans; in 
each scan the time on-source is of about 2.5 s. The high system 
temperature of the 8.3 GHz receiver was due to malfunctions in 
the receiver cryogenic system. 
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80 
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observed by the satellite, according to its publicly available 
pre-programmed pointing lis10. 

The SiMPlE project obtained 21 epochs of allocated 
time in the period between June-December 2010. Each epoch 
was scheduled for up to 24 hours. Only one epoch was com- 
pletely lost because of adverse weather conditions, which 
also affect 8.3 and 5 GHz observations. In the first six months 
of the SiMPlE project no 22-GHz acquisitions could be per- 
formed because of the insufficient sensitivity of the available 
K-band receiver. The temporary installation - for commis- 
sioning purposes - of a multi-feed 18-26 GHz receiver, whose 
final destination is the new Sardinia Radio Telescope, al- 
lowed us to carry out 22-GHz observations, together with 
new lower frequency ones, during the first semester of 2011. 

The scheduling procedure prioritised for each epoch the 
sources that were being observed by Planck, then the sources 
that in previous runs had been flagged out because of poor 
data quality, and, finally, all the other sources, in order to 
observe the sources as many times as possible. 

3.2 OTF cross scans and scheduling criteria at 
the Medicina radio telescope 

Observations were carried out with the 32-m single dish 
of the Medicina radio telescope in the OTF scan mode 
(Mangum et al. 2007), exploiting both hardware devices and 
software tools recently developed, at present still under com- 
missioning. On the hardware side, a new analogue backend 
was employed, the first one to be fully dedicated to contin- 
uum single-dish activities with the Medicina VLBI antenna. 
Besides permitting the execution of high-speed scans, is very 
helpful to avoid system instabilities and to better trace at- 
mospheric variations. Fast OTF scans were made possible 
as the antenna was provided with a new control system, the 
Enhanced Single-dish Control System (ESCS) specifically 
designed to perform single-dish observations exploiting the 
full potential of the telescope. Tests have demonstrated that 
scans can be carried out up to 20° /min without compro- 
mising the pointing accuracy. The optimised setup for the 
cross-scans required acquisitions at much lower speeds, with 
a sampling rate of 25 Hz and a resulting spatial resolution 
of 60 samples/beam. Table [1] lists the main scan setup pa- 
rameters employed. 

The number of cross-scans to be performed on each 
source was adjusted in real time during the observations. 

Informations on the Planck scanning strategy and 
pointing are publicly available for external observers at 
|http: //mmm. sciops ■ esa. Int/index .php?project=PLAHCK&page=Polntlng| 
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Table 2. Number of cross scans performed to achieve SNR=10 
on sources of the given fiux density limit, considering average 
weather conditions and the actual system setup exploited during 
the observing session. A minimum of 2 cross scans was carried out 
even on the brightest sources, despite a single scan would have 
been more than sufRciont from the sensitivity point of view, to 
cross-check for possible artefacts in the data. 



Flux density limit 
[m.Iyl 


5 GHz 


8.3 GHz 


200 


7 


22 


500 


2 


4 


1000 


1 


1 



according to the actual system performance and weather 
conditions, in order to reach a SNR of at least 10 in the 
final integrated scans. This dynamical scheduling was possi- 
ble thanks to an on-purpose-developed tool: the Positional 
On-the-flight-scan Planner (POP). As run-time execution 
of this tool completes in few seconds, POP allowed us to 
quickly produce new schedules whenever the weather and/or 
system conditions changed during an observing session. The 
only input required to the observer consists in the target co- 
ordinate list, the LST time at which the schedule is supposed 
to be executed, the observing frequency and, optionally, the 
estimated flux density of the sources. POP schedules the 
sources after checking for their visibility, taking into account 
the telescope movement limits and duration in azimuth and 
elevation, sorting them in ascending RA, to minimise the 
slewing time between targets and maximise the number of 
observed targets during the allocated time. 

For each source POP calculates the minimum number 
of scans required to reach a noise level 10 times smaller 
than the 23-GHz WMAP flux density, provided in the input 
position list. If no flux density is provided POP schedules 
for each target a number of scans equal to a user-deflned 
value indicated in the software configuration settings: this 
is of use in case of bad weather conditions or instrumental 
failures when a longer integration time is needed. 

Table[2]shows some examples of the cross scans required 
to obtain SNR=10 for different flux density limits. 

Each sub-scan on the sample sources is 5 Half Power 
Beam Width (HPBW) long, both at 5 and 8.3 GHz, in or- 
der to always have enough off-source samples and thus to 
better evaluate the background. In order to keep the sam- 
pling rate flxed at different frequencies (60 samples/beam), 
it required to change the scan speed, according to the dif- 
ferent HPBW. This means the Gaussian produced by the 
source is characterised by an high sampling rate, allowing 
to perform the fit on many more points. Scans on calibra- 
tors were set larger (7 HPBW): because of their crucial role 
in the sources flux densities recovery, the larger number of 
off-source samples provides a more precise identiflcation of 
the baseline, translating into a more reliable measure for 
the parameters which will subsequently be employed in the 
calibration phase (details in Section 

3.3 Region mapping 

Some sources in the NEWPS five-years 5ct catalogue are 
classified as undefined because of the lack of counterpart in 



low frequency catalogues. Some other objects have multiple 
associations within the beam size with low frequency wide 
area surveys (NVSS and/or GB6). These samples might in- 
clude spurious detections on WMAP maps as well as ex- 
tended sources or peaks of Galactic foreground emission er- 
roneously identified as sources by the detection procedures. 

In order to identify or to resolve the candidate sources 
we scanned patches of the sky around the source positions. 
We selected a sample of 29 undefined targets using coor- 
dinates of the most likely counterpart in the GB6 or NVSS 
catalogues when possible, otherwise we used the coordinates 
of the NEWPS detection. 

For each of these sources, we pointed and mapped patches in 
the sky of about 50' x 50', centred in the estimated position 
of the sources. The scanning strategy consists of 21 equally 
spaced scans in each direction (RA-DEC). Furthermore, con- 
secutive scans are performed with opposite movement direc- 
tion. To get a visual display of the mapped regions, we have 
built a grid based on the RA and DEC coordinates of the 
considered source and we associate the average value of a 
scan lying in each grid step to the corresponding RA, DEC 
position. The resolution of the grid on which the map is built 
could be changed, but in general, for a good visualisation of 
the sky patch, we decided for values ranging from 20 to 25 
grid steps. This is practically equivalent to a re-bin of the 
scans in steps ranging from 1.75 to 1.4 arcmin. 



4 DATA REDUCTION 

4.1 Data quality inspection and flagging 

In optimal conditions for observations, each scan has the 
shape of the beam transfer function, i.e. a Gaussian with a 
FWHM corresponding to the beam siz^H overlapped to a 
baseline that corresponds to the off-source zero level of the 
signal. Along the scans, the fluctuations of amplitude are 
given by a Gaussian noise. However, cloudy weather, the 
presence of random contributions by radio frequency inter- 
ference (RFI) or digital noise heavily affected portion of the 
data appearing as bumpy baselines or spike-like features. 

Hence, by fitting a Gaussian function and a linear base- 
line on the data it might be possible to assess the quality 
of each scan comparing the fitted value with the expected 
values and the overall goodness of the fit trough a anal- 
ysis. A similar test allows one to draw considerations about 
the Gaussian FWHM, the length of the baselines, the signal 
to noise ratio, the differences between the slopes of the two 
parts of the baseline (i.e. left and right with respect to the 
Gaussian). The broad variety of possible scan behaviours 
and the large number of parameters to be considered did 
not allowed to easily limit the parameter space. Further- 
more, in order to characterise the instrumental capabilities 
it is interesting, during the on-going commissioning phase, 
to classify the problems that affect the data. For this reasons 
we decided to manually inspect the whole dataset to classify 
the scans according to their quality and to the features that 
affect them, if any is present. Hence, we have developed a 

^ The measured deviations of the beam shape from a circular 
Gaussian beam are less than 2 per cent, including the measure 
errors. 
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Sample Number 



Figure 2. A sub-scan of a source (solid grey line) displayed by 
the flagging procedure. Here are also plotted the Gaussian fit of 
the scan (3 dot-dash line) and the linear fit of the left-right-hand 
baseline (dashed lines). 

tool that performs the fit of a Gaussian plus a linear base- 
line, showing it overlaid to the scan data as a reference (see 
Fig. [2)| in the visual inspection to guarantee homogeneity to 
the flagging criteria. 

We pay particular attention to the contributions affect- 
ing the linearity of the baseline: bumps, RFI, different slope 
of the left and right baseline arms. These features heavily 
affect the linear and the Gaussian fit of each scan and thus 
scans presenting any of these contributions are rejected. Of 
course any scan characterised by an odd shape of the beam 
transfer function is rejected as well. 

Bad weather conditions also reduce the SNR, so that in 
the worst cases (i.e. faint sources with high noise level) the 
sources were completely embedded in the noise fluctuations. 
Averaging over several scans reduces the noise and amplifies 
the SNR, thus for these critical sources the flagging was 
performed also on the integrated scans. 

In future phases we plan to automatise the flagging pro- 
cess on the basis of a statistical analysis of the flagging done 
so far, to minimise the subjectivity introduced in the flag- 
ging procedure by visual inspection. 

4.2 OTF scan processing 

We have developed the OTF Scan Calibration-Reduction 
(OSCaR) pipeline for the data reduction of calibrators and 
target sources. It constitutes an easily customisable ensem- 
ble of IDL routines capable of handling up to huge amounts 
of data, operating at all the SiMPlE frequencies. 

A run of the entire pipeline does not require a large 
amount of memory as well as a huge computational power: 
for about ten days of flagged data the complete process re- 
quires less than ten minutefl The entire software package 
can be easily modifled to be used for projects involving OTF 
cross-scan observations. 

^ On a machine equipped with a 2.0 GHz Intel Core 2 Duo and 
3 Gigabytes of RAM memory. 



The main steps of OSCaR (see Fig. [5]) are described in 
the following sections and can be summarised as follows: i) 
the estimation of the factor to convert the count sampling 
levels to Jy units by rescaling the counts measured on cali- 
brators to their known flux densities; ii) the reconstruction 
of the variation with time of the counts-to-Jy factor and of 
the component of the flux density error due to calibration; 
Hi) the calculation of the source flux densities and their er- 
rors, by applying the correct counts-to-Jy factor, integrating 
the scans for a source, and fitting them with a Gaussian that 
reproduces the receiver response function. 



4..2.1 Recovering of the counts to Jy conversion factor. 

The first task of OSCaR consists in recovering the conver- 
sion factors to transform the raw signal intensity, given in 
arbitrary counts, to a calibrated fiux density (Jy). This is 
done by fitting the scans with a Gaussian fit, plus a linear 
fit, according to the following: 

f{x)=Aoe^p{~z^/2)+A-i+AiX, (1) 
with 



In these equations, Aq is the height of the Gaussian, 
Al is the centre of the Gaussian, A2 is the width of the 
Gaussian, A^ is the constant normalisation term, and A4, 
is the slope of the baseline. In this way, we measure the 
Gaussian curve amplitude and compare it with the source 
absolute flux density. The fiux density/amplitude ratio is the 
counts-to-Jy factor valid for the observed elevation position. 

We consider as primary calibrators 3C286 and 3C295, 
the only two sources reporting negligible variability by Ott 
et al. (1993) with respect to the original scale of Baars et 
al. (1977). This flux density stability has been conflrmed by 
recent Effelsberg observations (Kraus & Bach, priv. comm.). 
Then, whenever it is possible, we measure the flux density 
of other calibrators (3C147, 3C48, 3C123, NGC7027, DR21) 
against the primary ones, exploiting the counts-to-Jy factors 
recovered on the primary calibrators. This allowed us the 
possibility of at least one calibrator observation per time 
interval of few hours (typically 2-3 hours in good observing 
conditions). 

A timeline of the conversion factors is recovered along 
all the days of acquisitions. In particular, once we collected 
all the conversion factors on an entire session of observation 
(see Fig.[3| we interpolated them in time, obtaining a contin- 
uous projection of these factors in the interval between two 
contiguous calibrator observations. We chose this strategy in 
order to have a time dependent conversion factor, reflecting 
in a smoother way the oscillations of each computed fac- 
tor during the day. These factors vary whenever a change 
happens in the weather conditions - implying a different at- 
mospheric absorption affecting the source ffux density - and 
when the system gain, for intrinsic reasons or user-defined 
choices, changes as well. Higher frequency observations are 
more sensitive to these contributions, and thus need a more 
frequent and accurate evaluation of the counts-to-Jy factors. 
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Figure 3. Conversion factor timeline recovered from an entire 
observing session. The cross and square symbols represents left 
and right channel respectively. 

4-2.2 Flux densities. 

The flux density measurements on the observed sources are 
performed through an algorithm which is very similar to the 
one used in recovering the counts-to-Jy conversion factors. 
The equation for the Gaussian fit is also the same. The ac- 
cepted scans are selected and integrated according to the 
tags provided during the flagging phase. The integration is 
performed through a precise positional alignment, in turn 
allowed by the very reliable antenna pointing system. It is 
also possible to define the time interval within which more 
scans on the same source must be considered consecutive and 
thus integrated, otherwise they are considered as separated 
measurements, each processed with the proper calibration 
factor. 

Before proceeding with the flux density computation a 
further check on the HPBW of the scan is performed. We 
accepted scans with values of the HPBW within 20 per cent 
variation from the nominal value. This criterion selects not 
only the scans for which the instrument behaved properly, 
but can also remove from our analysis sources which are ex- 
tended by more than 1.2 x HPBW, for which our flux density 
estimates are only a lower limit. 

We obtain separate flux density estimates for the two 
acquisition channels and for each scan direction. We firstly 
compute the flux density from the cross created by the inte- 
gration of the RA-DEC scans in each channel, then through 
a weighted average we achieve a unique flux density for the 
considered source. The application of the counts-to-Jy fac- 
tor takes into account that the source and the calibrator had 
been observed at different elevations. The conversion factor 
is properly rescaled by means of the standard gain-elevation 
curve provided for the Medicina dish. The equations are the 
following: 

G = (-4.6834953 ■ 10"^ ■ e^) + (3) 

-(-(6.2403816 ■ 10"^ ■ e) + 7.9212981 • 10"^ , 
at 5 GHz, while at 8.3 GHz it becomes 
G = (-7.2457279 ■ 10"^ ■ e^) + (4) 

-^(1.0623634 ■ 10"^ • e) + 6.1059261 • 10"^ . 



,150 

-150 -100 -50 50 100 150 

A RA [arcsec] 

Figure 4. The difference between the input and the recovered RA 
and DEC position of the Northern Sample sources for about 200 
individual pointings at 5GHz (for comparison the beam FWHM 
is 7.5 arcmin). 

G represents the elevation dependent normalised gain, while 
e is the elevation of the considered source/calibrator. A 
flux density compensation related to a possible offset in the 
pointing of the sources is also taken into account according 
to the following equation: 

= exp[~{x ■ 1.66/ {HPBW ■ tt))] (5) 

In this equation j — RA, DEC, k = L, R indicating left 
and right channel, x represents the positional offset accord- 
ing to the direction considered, and HPBW is the nominal 
value of the beam at the frequency considered. We calcu- 
late the positional offset by comparing the position of the 
peak detection performed by the fitting function with the 
input coordinates values. The positional offset of about 200 
individual pointings is shown in Fig. [4] 

Finally, each source is calibrated using a conversion 
factor recovered matching the acquisition epoch to the 
calibration timeline. In this way we ensure a short time 
distance between the observed source and the calibration 
factor, in order to minimise the possibility of weather or 
system variations. 

The set-up of the run is interactively performed after 
launching the main program. Among the available options, 
there is the possibility to call each procedure independently, 
to provide a list of sources to be analysed or to work on an 
entire observing session. 

Through the tags provided during the flagging phase we 
analyse each channel of each sub-scan independently. While 
this procedure would be standard for a polarisation analy- 
sis, for the total intensity study we have the possibility to 
choose to use one or both channels of each sub-scan. 
Particular attention was dedicated to the handling of the 
time tag of each scan. We aimed at collecting as many obser- 
vations as possible, of each source, at different times. These 
time intervals can span from hours to weeks or months. Thus 
we implemented the pipeline in order to have a coherent in- 
tegration of the scans concatenated in time. The automati- 
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Figure 5. Scheme of the main operations in OSCaR. See 14.21 for 
a detailed description of the procedures. 



sation of this strategy of integration permits to easily handle 
multiple observations of calibrators in a single day and con- 
sequently the creation of the detailed calibration timeline 
we use to calculate the flux densities. The possibility to re- 
cover the flux densities of selected sources in more observing 
sessions at once makes OSCaR well suited to evaluate vari- 
ability as well. 

We could have relied on the telescope point- 
ing/acquisition system for an immediate, sample- by- sample 
integration of the scans. Through a dedicated test, we con- 
firmed that the maximum misalignment encountered, com- 
paring the initial positions of the scans, corresponds to one 
sample of data (7.2 arcsec for 3°/min scans, i.e. less than 2 
per cent of the HPBW). However, for a more precise estima- 
tion of counts-to-Jy conversion factors and flux densities, it 
is important to avoid any possible offset, thus any misalign- 
ment, even of such entity, is corrected during the integration 
of the scans. 



4-2.3 Flux density error estimation 

The error budget in the final flux density results from two 
main components: the calibration error Seal and the noise 
error Sn 



The two terms have identical physical causes, although 
with different weights, and can be treated using a common 
formalism. The contribution to each component is given by: 



~ V ^baseline + ^Gauss 



(7) 



(6) 



with i = cal, N. 

The basehne error (SbaseUne) is mainly due to the white 
instrumental and background noise and to the varying in- 
fluence of the ground radiation while the antenna is scan- 
ning at different angles, ft has two spectral components: one 
showing white spectrum (background, atmosphere, statistic 
uncertainty) and the other characterised by 1/f behaviour 
(long time drift of detectors and slow atmospheric changes 
between the telescope and the target). The confusion limit, 
the asymptotic limit of the white component in the case of 
a large number of samples, represents only a part of the 
baseline error. 

In order to compute this contribution we subtract the 
linear fit, performed on the scan neglecting the Gaussian 
bell, from the baseline. This reduces the long time drift ef- 
fect. Then we calculate the standard deviation of the sub- 
tracted data. In this way, we take into account both the 
white noise contribution and the spurious residual effects on 
time scale of a single-source observation. As a general rule, 
we decided to keep margins, when fitting the linear base- 
lines, neglecting the first ten points on the left-hand and 
on the right-hand sides of the Gauss fit, whose width was 
calculated from the beam size. 

The Gauss error {Sgauss) estimates the goodness of the 
Gaussian fit. It represents the error committed in calculating 
the coefficient (amplitude of the Gaussian) of the expo- 
nential term, and can be ascribed either to white scatter 
effects affecting the whole scan or to local drops and spikes 
of signal, mostly in correspondence of the tip. However, in 
most cases, the preliminary operation of data flagging is ex- 
pected to have already prevented the Gaussian shape from 
evident anomalies. Operatively, the Gauss fit is calculated 
over the integrated scans and the associated error is pro- 
vided by the fitting routine itself. The final error in the fiux 
density, indicated in Eq. ((7|), comes out as the amplitude 
uncertainty of the height of the bell estimated through the 
Gauss fit. 

Because of the flagging strategy adopted, acting individ- 
ually on each subscan, the number of good data was different 
for RA and DEC scans as well as for the two polarisation 
channels, for each source observed. For this reason, Eq. (0 
was calculated also for right-hand channel (RHC) and left- 
hand channel (LHC) when scanning in RA and in DEC. 

In accordance with the observing strategy, about 20 
sources were observed in the time window between two cal- 
ibrator^. In a few cases the primary calibrators were not 
available during the whole observing session: consequently, 
we have been forced to calculate the calibration constants 
using secondary calibrators. Through dedicated tests, we 
found out that this strategy could affect the calibration er- 
ror by less than 5 per cent in a day with good observing 
conditions. Therefore, we decided to calculate a weighted 



^ Performing an OTP cross-scan on a source takes approxima- 
tively 30 sec, including pre and post scan operations. 
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calibration error for each source, relating the calibration er- 
ror to the time window considered instead of calculating an 
average error over the whole run. 

The calibration error on each observed source comes out 
from weighting the errors committed on the two nearest cal- 
ibrators by the distance in time of each source from the 
calibrators. Hence, given ri and T2 the distance in time of 
the i-th source from calibrators 1 and 2 (characterised re- 
spectively by errors AScai^^nd AS'caia)) ^^e weighted error 
is: 



AScal-^ , A,Scaln 
2--CAL ~ T2 ; h Tl - 



Tl + T2 



Tl + T2 



(8) 



Similar to ASi, T,cal was calculated for channels RIfC 
and LHC both in RA and DEC; when combining these error 
terms with Eq. © we get the following four-term array: 



ASra, 



ASi, 



E2 

CAL 



RA, + ^SIr^ 



%AL' 



ASdeCl ~ 



ASdeCr, = 



5^2 

CAL^'^^L 



+ AS'2 „ „„ 



CAL 



DEC„ + ^SlnECr, 



(9) 
(10) 

(11) 
(12) 



The calibration errors, calculated following the above 
procedure, were hence compared with the error calcu- 
lated just considering the whole error bar obtained when 
displaying the calibration constants referring to all the 
calibrators available during the observing day considered. 
This comparison is displayed in Fig. [S] 

The final error T.fx in the flux density of each source comes 
out as the square root of the sum of the squares (RSS). 



(13) 



The error contribution can be mapped through a plot, 
showing how it varies (in percentage) with the source inten- 
sity (Fig. [3, at 5 and 8.3 GHz. 



4.3 Simultaneous ATCA-Medicina observations 

In addition to the NEWPS Northern hemisphere sample, we 
have observed with the Medicina radio telescope a list of 12 
sources with 5 < 0° among the targets of the PACO project. 
We observed almost all these targets simultaneously with the 
two facilities. In any case, the observations were performed 
within 10 days from the Planck satellite ones. The PACO 
project observed them in 24x512 MHz-wide frequency sub- 
bands between 4.5 and 40 GHz, including data at 5.244 and 
8.232 GHz. This sample, although small, has been used as 
a test to verify that our procedures give consistent results. 
The two projects in fact are characterised by the use of dif- 
ferent telescopes, data-reduction pipelines, and completely 
different calibration schemes and sources. 

The SEDs of the equatorial targets including the Medic- 
ina and ATCA radio telescopes, and the Planck satellite data 
are in Fig. [§1 showing a fairly good agreement. 

The comparison is unbiased by source variability and 
only depends on the instrumental properties. The Planck 
satellite lower frequency channel is too high to transfer 
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Figure 6. A typical conversion factor timeline. Tiie black solid 
line represent tiie effective conversion factor applied to the data, 
while the grey zone represents the error associated with the fac- 
tors themselves. The dashed lines represents a rough estimate of 
the error that can be considered for the whole day of year (DOY). 
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Figure 7. The flux densities from 100 measurements of the 
Northern Sample sources (crosses) versus their percentage error 
are plotted. The upper panel shows the 5-GHz measurements, 
while the lower those at 8.3 GHz. The black lines represent the 
fits of the plotted data and the same quantity for the other fre- 
quency is overplotted for comparison (dashed lines) . Multiple ob- 
servations of the same sources are also considered. 
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Figure 8. SEDs for a sample of equatorial targets observed simultaneously with the Medicina radio telescope (filled diamonds), the ATCA 
(cross, in the framework of the PACO project) and the Planck satellite (empty diamonds, data from ERCSC, The Planck Collaboration 
2011a). The target 0238-084 was observed from Medicina with some days of delay. The displacement in its SED could be due to variability 
of the source itself. 



its CMB dipole based flux density calibration to the low- 
frequency Medicina telescope channels, but we can compare 
the two ground-based facilities flux density scales. The flux 
density calibration for the ATCA data is based on the obser- 
vation of one single very stable source, PKS B1934-638. Fig. 
[5]shows the flux density comparisons at 5 and 8.3 GHz. The 
best fitting line at 5 GHz has a slope equal to 1.03 ± 0.09 
and crosses the y-axis at (—0.39 ± 0.36) Jy. At 8 GHz it 
has slope equal to 1.02 ± 0.15 and crosses the y-axis at 
(—0.62 ± 0.65) Jy: the flux density is consistent within the 
error bars between the two instruments. A small bias in cal- 
ibration is probably present, but the small sample does not 
allow to quantify it. Further observations of a larger equato- 
rial sources sample are on-going to quantify the calibration 
differences between the telescopes with smaller uncertain- 
ties and to have a common calibration scale between several 
facilities and including at mm frequencies also the Planck 
satellite. 



5 THE 5 ^ 45° SAMPLE 

After the flagging and the removal of data points with signal- 
to-noise ratio lower than 5, 67 of the 69 NEWPS selected 
sources with 5 45° observed with OTF cross-scans have 
reUable data at 5 GHz (the exceptions are the 179th and 
the 245th in the NEWPS catalogue for which all the data 
were flagged for bad weather conditions) and 59 sources have 
valid data also at 8 GHz in at least one epoch in the June- 
December runs. Table 3] lists the data at 5 and 8.3 GHz for 
one observing epoch for each source per row. The columns 
are as follows: 




Figure 9. Flux density comparison between SiMPlE and PACO 
data at 5 and 5.244 (filled symbols) and at 8.3 and 8.232 GHz 
(empty symbols) respectively. No corrections for the small differ- 
ence in frequency have been applied. The best fitting lines at 5 
and 8 GHz are respectively the solid and dashed lines. The dotted 
line represents unity. 



1 Source identification. For sources in the NEWPS cata- 
logue the id corresponds to the sequential number in the 
Massardi et al. (2009) catalogue. 

2-3 Right ascension in hour and declination in degrees. 

4 Epoch of observation in YYYY-MM-DD. 

5 Flag asterisks for simultaneity with the Planck observa- 
tions. 

6-7 Flux density at 5 and 8.3 GHz in mjy. 
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8-9 Flux density error at 5 and 8.3 GHz in mjy. 

10 more targets have been observed by mapping the 
surrounding sky region and will be discussed in the following 
section. 



5.1 Region maps with 5 ^ 45° 

Among the regions that we have mapped at 5 GHz, ten have 
declination S ^ 45° . The maps are shown in Fig. 1111 

For each map we perform a two-dimensional Gaussian 
fit, in order to identify the candidate source and estimate a 
lower limit to the source flux densities. This procedure fits 
an elliptical Gaussian equation to gridded data. The fitting 
function is: 

F{x, y)^Ao + Ai- exp(-t//2) , (14) 

with the elliptical function described by ?7 = (x'/a)^ + 
(y'/b)'^, while Ao and Ai are the constant term and the 
scale factor, respectively. The rotated coordinate system is 
defined as: 

x' = - /i)cos(r) - - fc)sin(T) (15) 
y' ^ {x - h) sm{T) + {y - k) cos{T) (16) 

In our case, x and y are the scan directions RA and DEC, 
{h, k) is the centre of the ellipse, and T is the rotation angle 
from the x axis in clockwise direction. Furthermore, it is 
possible to recover the width of the Gaussian in the two 
directions, the location of the centre on the two axes, and 
the rotation of the ellipse with respect to the x axis. Once 
the peak position has been identified we estimated the peak 
flux density which is a lower limit to the source flux density 
in case a candidate could be identified. As the profile of an 
extended source could be far from a Gaussian this approach 
has been used only to identify a peak and estimate a peak 
flux density. As the peak flux density is only a lower limit 
of the integrated flux density, we define a candidate source 
detection only the case in which the peak flux density is at 
least 3 times larger than the noise estimated on the map and 
the source can be univocally identified in the map (i.e. there 
is only one peak and the scans passing across it passed the 
visual flagging procedure). Table[3]lists the fit details for the 
10 maps in the Northern sample region. 

5 of the 10 mapped regions with S ^ 45° have been 
identified as candidate sources, but only three of them have 
(peak) signal-to-noise ratio larger than 5. Fig. [10] shows a 
graphical elaboration of the map of n22, one of the five 
source detections. This object has been identified with the 
HH region NGC 0281. The recovered peak position of the 5 
detected sources was then compared with that provided by 
WMAP: in all the cases the displacement does not exceed 
10 arcmin. 

Hence among the 79 sources identified in the NEWPS 
sample with S ^ 45°, about 95 per cent have been con- 
firmed as genuine/candidate sources. A proper analysis of 
the comparison with the selection flux densities will be done 
in future works where we will present the 22-GHz data for 
our objects, for which the observing campaign is on-going. 



Table 3. Peak flux densities extracted for the 10 mapped regions 
with 5 ^ 45°. The last column indicates as 1 the sources that 
have Speafc SG-ffz/"" > 3 and the scans across the peak position 
passed the visual flagging procedure. 
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Figure 10. Surface projection of the sky region centred on 
NEWPS source n22 identified with the HII region NGC 0281. 
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Figure 12. Distribution of spectral indices in the 5-8 GHz fre- 
quency ranges. 
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Figure 11. The ten regions of sky mapped around the undefined sources described in l3.3l The flux density scale is expressed in Jy. 



The Simultaneous Medicina- Planck Experiment 13 



20 




Figure 13. Distribution of spectral indices in the 5-30 (shaded 
area) and 30-70 GHz (hatched area) frequency ranges. The bin 
size is fixed at 0.5 for both the histograms. The displacement is 
due to the different ranges characterising the two spectral indexes. 




Sgb6 [Jy] 

Figure 14. Comparison of SiMPlE and GB6 5-GHz flux densi- 
ties. The solid line shows the best linear fit and the dot dashed 
lines enclose the area between the best fit plus and minus one 
standard deviation of the flux density ratio. 



5.2 Spectral behaviour and comparisons with 
other catalogues 

For each source we have identified the best epoch of obser- 
vations as tlie one with both 5 and 8.3 GHz observations. 
In case of multiple epochs with both frequencies we have 
chosen those having higher signal-to-noise ratio at both fre- 
quencies. We thus have a sample of 61 sources for which we 
estimated the spectral indices af "^ (defined according to the 
flux density being cx i^"). The median spectral index is -0.09 
with a standard deviation of the distribution equal to 0.62 
(see Fig. I12p . This is compatible with the 5-10 GHz spectral 
index found for the PACO bright sample (the selection was 
based on 20-GHz flux densities > 500 mjy). 



33 sources have counterparts within the Planck ERCSC 
(The Planck Collaboration 2011a). The median of the 5-30 
GHz spectral index is -0.07 (see Fig. ll3|l . As elsewhere stated 
(Massardi et al. 2008, 2009, 2011) the bright high-frequency 
selected samples are mostly constituted by flat spectrum 
sources in this frequency range (see also The Planck Collab- 
oration 2011c, d). 

For this sample in the ERCSC the 30-70 GHz spectral 
index is -0.16. As compared with the 100-GHz flux density 
limited sample in Sadler et al. (2008) , the steepening is more 
pronounced for brighter samples. 

We have compared our sample flux densities with the 
4.85 GHz GB6 catalogue (Gregory et al. 2006). 60 objects 
in the Northern sample have a counterpart in the GB6 cata- 
logue (the remaining objects include the 10 sources for which 
we have performed maps of the surrounding region and the 
objects in the region with 5 > 75° that is not covered by 
the GB6 survey). The flux density comparison is shown in 
Fig. 1141 Performing a linear flt of these data we found the 
following relation between the flux densities from SiMPlE 
and GB6: 

Ss^MPlE,5/[3y] = (0.900±0.057)-Sgs6 + (0.250±0.130).(17) 

For this wide sub-sample the rms of the fractional dis- 
placement is 43 per cent: since more than 14 years divide the 
two catalogues, this effect is probably entirely due to source 
variability. This is comparable with the 38 per cent variabil- 
ity found, at the same frequency on few years timescales, 
on the PACO sample confirming earlier indications that the 
variability increases with the time lag, for lags of up to sev- 
eral years. Our results seem to indicate that the asymptoti- 
cal value for the source variability at 5 GHz is ~ 43percent. 



6 SUMMARY 

In this work, we have described the sample selection, the 
observing strategy, the data reduction procedures, and the 
software tools developed for the Simultaneous Medicina 
Planck Experiment (SiMPlE) project. The selected sample 
includes 253 sources in the NEWPS sample and 11 candi- 
date variable Galactic sources. 

We have presented some tools that we have developed 
to analyse data collected with the Medicina 32-m single dish. 
The tools allow us to efficiently schedule source lists, to visu- 
ally inspect the scan quality, to compose and analyse maps 
of sky patches, and to calibrate and estimate the source flux 
densities. 

The developed pipeline allowed us a fast and detailed 
automated analysis of our entire data-set, and the whole 
software package can be easily adapted for analogous OTF 
scan-based experiments. 

As a result we have shown the results obtained for a 
sub-sample of 79 sources with S ^ 45° observed in 22 runs 
in the June-December 2010 epoch, at 5 and 8.3 GHz. Data 
for 2 of them have been flagged for bad weather conditions. 
For 10 targets we have performed maps of the surrounding 
sky area conflrming the detection in the WMAP maps for 4 
of them. These results conflrm that the NEWPS catalogue 
could be considered 95 per cent reliable. 

The comparison with the ERCSC and the analysis of 
the spectral behaviour have conflrmed that bright sample 
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sources are mostly composed by flat-spectrum objects up 
to 30 GHz, but show an overafl steepening in the spectrum 
in the frequency range 30-70 GHz. The comparison of the 
60 sources in common with the GB6 catalogue sliow a rms 
fractional displacement due to source variability of 43 per 
cent that is an indication of the long-term variability for 
bright sources on time scales of more than 10 years. 

In addition to the SiMPlE Northern Sample, we also 
report the observations for a list of 10 sources with 5 < 0° 
chosen from the targets of the PACO project. The compar- 
ison of the flux densities obtained with the two instruments 
and the SEDs of these targets showed that the ATCA and 
the Medicina radio telescopes have a consistent flux den- 
sity scale, despite the fact that the calibration and the data 
reduction follow independent methods. Larger samples are 
being observed simultaneously with the two telescopes to 
confirm this finding. 

Observations for the SiMPlE experiment are continuing 
for the first semester of 2011 in order to overlap with at 
least 2 full surveys of the Planck satellite and exploiting also 
the commissioning phase of the 22-GHz multi-feed receiver 
at the Medicina radio telescope. The SiMPlE sample might 
constitute a reference sample for this new telescope and for 
on-going Northern Hemisphere surveys. 
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Table 4. A sub-set of the SiMPlE Northern sample. The table will be replaced by the complete sample. 
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